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P\[ , Abstract. This paper is concerned with the traveling waves of de- 

layed reaction-diffusion systems where the reaction function possesses the 

^ [ mixed quasimonotonicity property. By the so-called monotone iteration 

scheme and Schauder's fixed point theorem, it is shown that if the system 
has a pair of coupled upper and lower solutions, then there exists at least 
a traveling wave solution. More precisely, we reduce the existence of trav- 
eling waves to the existence of an admissible pair of coupled quasi-upper 
and quasi-lower solutions which are easy to construct in practice. 
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1 Introduction 

Reaction diffusion system is used to model the spatial-temporal pattern. In 
the past decades, the traveling wave solutions of the reaction diffusion systems, 
which are studied as a paradigm for behavior, have been widely investigated due 
to significant applications in chemical engineering, population dynamics and 
biological models. Since the first instances in which traveling wave solutions 
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were investigated were given in 1937 by Kolmogorov et al. [H] and Fisher p], 
many methods have been used to study the travehng wave solutions of various 
parabohc equations and systems, for example, the phase plane technique in 
IIIEIIIIIIHEZIEIESIES], the degree theory method and the conley index 
method developed in P E El [S] ■ 

In many realistic models, the delays should be incorporated into the reaction 
diffusion system. Due to the presence of delays in the reaction diffusion system, 
the classical phase plane technique can not generate a monotone flow which 
ensures the existence of the traveling wave solution. Recently the classical 
monotone iteration technique was first used by Wu and Zou [361 EI] to establish 
the existence of traveling wave solution for delayed reaction diffusion system. 
They employed the idea of upper and lower solutions and an iteration scheme to 
construct a monotone sequence of upper solutions which was proved to converge 
to a solution of the corresponding wave equation of the reaction-diffusion system 
under consideration (see also [9l HO]). In fact, many researchers had used the 
monotone iteration technique to prove the existence of the reaction diffusion 
system in [HJlSl [ISl [ISl [HI [IHESlEHl EH]. In [231 El], Ma et al. proved 
some existence results for traveling wavefronts of reaction-diffusion systems by 
using Schauder's fixed point theorem. One important feature of Ma's method, 
which was different from the work of Zou and Wu [36], was that the upper 
solution of the wave equation was not necessary to converge to two distinct 
trivial solutions when t — )■ — oo and t — > +cxd respectively. Li et al. developed a 
new cross iteration scheme and established the existence of the traveling wave 
for Lotka- Volterra competition system with delays [191 EHl EI] • 

More recently, Boumenir and Nguyen discussed in P] a modified version of 
Perron Theorem for C^-solutions, and set up a rigorous framework for the mono- 
tone iteration method and then apply it to the predator-prey and Belousov- 
Zhabotinskii models with delays. 

However, in the iteration process by the monotone iteration it is required 
that the nonlinear reaction function possesses a quasimonotone property in the 
sector between the upper and lower solutions [21 ESI |36] . This paper focus on 
the delayed reaction diffusion system without quasimonotonicity. Motivated 
by the above work and the upper and lower solution method developed by Pao 
[SHlEHlEOlET] and Li et al. [19] for reaction diffusion systems, we use the coupled 
upper and lower solutions to deal with the non-quasimonotonicity, which was 
first given out in [22j. Via the coupled upper and lower solutions, we construct 
an appropriate closed bounded convex set. By use of the Schauder 's fixed point 



theorem in the convex set, we show the existence of the travehng wave solution. 
Moreover we reduce the existence of travehng wave solution to the existence of 
an admissible pair of quasi-upper and quasi-lower solutions which are easy to 
construct in practice. 

This paper is organized as follows. In Section 2, we show the existence of the 
traveling wave solution by constructing the classical coupled upper and lower 
solutions. In Section 3 we relax the classical coupled upper and lower solutions 
to the C^ smooth coupled quasi-upper and quasi-lower solutions. Section 4 deals 
with systems with quasimonotone nondecreasing functions, and the definition 
of ordered quasi-upper and quasi-lower solutions is introduced and an existence 
result of a traveling wavefront is given by the monotone iteration method. In 
Section 5 the main result is illustrated by and applied to a delayed Belousov- 
Zhabotinskii equation and a Mutualistic Lotka-Volterra model. This paper ends 
with a short discussion. 

2 Coupled upper and lower solutions 

In this paper, we will consider the following system of react ion- diffusion systems 
with time delays 

d d^ 

—u{x,t)-D—u{x,t) = f{u,Ur), (2.1) 

where x G M, t G (0, oo), u = (ui, ■ ■ ■ , u„) G M", u^ = (ui(x, t—Ti), ■ ■ ■ , Un{x, t— 
Tn)) G M" for some positive constants n, ■ ■ ■ ,r„, which are so-called discrete 
delays. and D = diag{di, ■ ■ ■ ,dn) with rfj > 0, f : R" x R" — )■ R" is Lipschitz 
continuous. 

For convenience, we denote by Cft(R, R) the space of all bounded and con- 
tinuous functions /i : R — )■ R endowed with the super-norm. Moreover, for any 
k G R+, we denote by C^(R, R) the space of all continuous different iable up 
to the [fc]-order functions h such that d^h G Cb(M, M) for any I7I < [k] {[k] 
denoting the integer part of k). The above spaces for vector- valued functions 
(with n components) are denoted by C5(R, R") and C^(M, R"), respectively. 

A traveling wave solution of (12.11) is a special translation invariant solution 
of the form u(x, t) = (p{t + x/c), where Lp G C^(R, R") is the profile of the wave 
and c > is a constant corresponding to the wave speed. The vibration at the 
space point a; = is u(t) = v^(t); The vibration u(t) propagating from the space 
value a; = to X costs the time x/c, where c is wave velocity, in the case c > 



the traveling wave move to the left, in the case c < the traveling wave move 
to the right. Our definition is more visual than u(x, t) = (p{x + ct) in [21 [23tl36]. 
Substituting u(x, t) = ip(t + x/c) into (12. ip and letting s = t + x/c, denoting 
also t, we obtain the 

ip'it) - ^ip"{t) = %(t), ^{t - r)), t G R. (2.2) 

corresponding wave equations If for some wave velocity c, (12. 2 p has a solution 
if defined on M such that 

lim (f{t) = u_, lim (f{t) = u+ (2.3) 

t— >■— oo t— >+oo 

exist, then u(x, t) = (f{t + x/ct) is called traveling wave with speed c. Moreover, 
if (fi is monotone in t e M, then it is called a traveling wavefront. 

Without loss of generality, we can assume u_ = and u_|_ = K > 0. Let 

C[o,K](M,M") = {^ = (<^i(t),--- ,<^„(t))Ga(M,M"): < v?,(t) < K,,teR}. 

Our aim is looking for a solution of (12. 2p in C[o,k](R;R")- Throughout this 
paper, the following hypothesis will be imposed on the reaction term f : 

(ifi) f (0) = f (K) = 0. 

Obviously, we should replace (12. 3p with 

lim ip{t) = 0, lim (p{t) = K (2.4) 

t—^—oo t—^+oo 

In this paper, we explore the existence of the traveling wave solutions of (12. ip 
where the reaction term f is mixed quasimonotone. 

{H2) The function f (u, u^) = (/i(u, u^), ■ ■ ■ ,/„(u, u,-)) is a C^ function and 
possesses a mixed quasimonotone property in a subset [0, K] of M". 

The above hypothesis implies that there exist constants /3j such that fi satisfies 
the Lipschity condition 

|/i(u,u^) -/i(v,v^)| < /3i(||u- v|| + ||u^ - v^ll) 

for all u, V, u^ and v,- in C[o,k](R, M"), z = 1, ■ ■ ■ ,n, where | ■ | and 1 1 ■ 1 1 denote 
the super norm in R" and C(R, R"), respectively. 



Recall that by writing vectors u and u^ in M" in the split form 

where [u]o- denotes a vector with a components of u, the function f (u, u^) is 
said to have a mixed quasimonotone property if for each i = 1, ■ ■ ■ ,n, there 
exist nonnegative integers at, bt, Ci and di with a, + 6j = n — 1 and Ci + di = n 
such that /i(u, v) is monotone nondecreasing in [u]a- and [vj^ and is monotone 
nonincreasing in [u]h- and [v]^.. li bi = di = for all i then f (u, v) is said to be 
quasimonotone nondecreasing. 

The above general assumptions are used to establish the existence of a trav- 
eling wave solution to (12. 2p . Our approach to the problem is by the method of 
coupled upper and lower solutions which are defined as follows: 

Definition 2.1. A pair of vectors ip = [ipi,--- ,ipn), = {0i,''' ,0n) in 
C^(M, R") are called coupled upper and lower solutions of Ii2.^) if (p > and if 

0'iii) - %0'lit) > fii.0h [^]a,:, [0\b,, [(prV,, [^r]dj, ^2 5) 

0'i{t) - %0'l{t) < fi{0U [0U, [0]b,, [0r\c,, [VrW (^ = 1, ■ ■ ■ ,n), 

where (prit) = ip{t — r). 

Remark 2.1. Iff{ip,ipr) is quasimonotone nondecreasing, that is, bi = di = 
for all i, then the pair of vectors called ordered upper and lower solutions of 

Since that f satisfies (-^^2) and Lipschitz continuous, we have 



fA^i^ Ma,, Mb,, Wr\c,, WM) - fA^'i. Ma,> Mb.' br]c,, Wr\d, 

+ P>i{^i - v^-) > for all < (y9- < (y9j < ii'j, i = 1, ■ ■ ■ , ra. 



(2.6) 



Next we define an operator H : C[o,k](K, M") -^ C[o,k](K, M") by 

H((^,(^,)(t) = f((^,(^,) +/3(^(t), ^ e qo,K](M,M"), (2.7) 

where H = {Hi, ■ ■ ■ ,if„), Hi{ip,ipr){t) = fi{(p,(pr) + I3iipi{t). Clearly, with the 
above notations, the system (12. 2 p is equivalent to the following system of ordi- 
nary differential equations 

ip'{t)-^ip"{t) + f3^{t)=U{^,^rm, tGM. (2.8) 



Our first iteration involves the following linear system of ordinary differential 
equations 



c(^(i)v_*^-(i)v'^ «.-(!) 



7 - f (^rO" + m' = P^^^ + /.(^., [^]a„ [^l^,, [^.]c., [^.]dj, ,^ g^ 



Note that 



c2(i-0T4M7^) , _c2(i + v/rT4M7^) 



are the negative and positive real roots of the equation 

^A2-A-A = 0, z = l,2,---,n. 
Using the Perron Theorem yields 

^^^^ = TnT^-XT^ r e^"(*~^)(A^, + /,(^„ [^],,, [^],,, [^Ac, [0rU))ds 

O-iyMi — Mi) J-oD 

/ + 00 
e^'^^'-'\Pm + M^i, [^]a,,, [0]b,, [^rU, [0rU))ds), 

. (2.10) 

^^ = Tn r^( / e^"^*"'^^^. + /.(^., [0U, i^k, [0rU, [^rU))ds 

"i(,/\2i — Mi) J-oo 

e^^'(*-^)(A<^i + fii0i, [0U, [<fk, [0rU, [<frU))ds) 

for z = 1, ■ ■ ■ , ra. Then by Lemma 2.1 of [22], x'-"'^^ = {x[ ,■■ ■ ,x^n^) and x*-"*^-* = 
{x\ \- ■ ■ ,Xn) have the following properties. 

Lemma 2.1. Let x^^' andi6^'> he the solution of \2. 9\) . then we have 

(l) (p < x(^) < x(^) < If. 

(ii) x*-^-*, x''"'^^ 0,''^^ 0. pair of coupled upper and lower solutions of Ii2.^) . 

Now in order to prove the existence of the traveling wave solution, we are 
in the position to apply the Schauder's fixed point theorem. We define the 
operator 

F = (Fi, ■ ■ ■ , i^„) : C[o,K](M, M") -^ C[o,K](ffi, »") 



by 



for i = 1, 2, ■ ■ ■ ,n. 



Let p > be such that p < min{— Aij, A2i : i = 1,2,- ■ ■ , n}, and let 

Bp(M,M") = {^e C[o,K](K,K") : suptm\f{t)\e-P^'^ < oo}, 
|(^|p = supigK|(^(t)|e-''l*l. 

Then it is easy to check that Bp(]R, M"), | ■ |p is a Banach space. 

Lemma 2.2. Let the closed convex set F = {</) G C[o,k](IR, I^") : <^ < V^ < <^}, 
where ip and are coupled upper and lower solutions of li2.S\) . then F{T) C T. 

Proof. Since that the mixed quasimonotone property of the operator H, \/ip G F, 



we have 

{F^^^m < T7V^^^( [ e^"^*"'^^^. + M'fi, [^U [0k, [0rU, [0rU))di 



oo 

+ 00 



+ / e^'^^'-'\/3,0i + fii^i, [0U, [0]b,, [0rU, [0r]dj)ds) 



By the virtue of Lemma [2.11 we induce {Fiipi){t) < (pi{t). Similarly, we have 
{Fi^i){t) > 0i{t) for alH = 1, ■ ■ ■ , n. Therefore F(F) C F. D 

Lemma 2.3. If the hypothesis (ifs) holds, then F : C[o,k](K,IR") -> C[o,k] (K, K") 
is continuous with respect to the norm \ ■ \p in Bp(R, R"). 

Proof. Vy^, ip' G C[o,K](IR5lR"')! in view of the definition of F, we have 



C2 



ai[Mi — Mi) J-oo /2 ii\ 

/+00 *- ' ' 

e^-(*-^)(A(v.. - if[) + f.{cp,cpr) - f.iip',if',))d.s). 

Since that ip, ip' < K, {H2) implies that /i(v?, fr) is bounded for (p G C[o,k](I^, I^" 
Then the term (3i{(pi — (p[) + /^(v?, V'r) — /ilv^', V'r) is bounded, for convenience, 
we denote [3i{ipi — (p[) + /i(v5, v^r) — /i(v^', v't) — ^i- (12.111) is transformed into 



^2 nj, rt /.+00 

c^Mi 1 1 c^M- 



(2.12) 



di{X2i — Aij) A2J Aij <iiAiiA2i 



It follows from (KT^ that 



F,(^,^.) -F,(^',^:)|e-^l*l < -;T^e-''l*l < --^. (2-13) 



diXliX2i Cl'iMi^2i 

Therefore \Fi{(p,(pr) — Fi{(p',(p'^)\p < — ^\^^{ for all i = 1, ■ ■ ■ ,n. That is, 
F : C[o,K] (M, K") -^ C[o,K] (K, K") is continuous'' D 

Lemma 2.4. // i/ie hypothesis {H2) holds andF is defined in Lemma lKE, then 
F : r — !■ r zs compact. 

Proof. First we compute |^|p, for any cp E C[o^k](^,'^"'), we have 

2\ 2t It J J — 00 / r\ 1 /I \ 

0'i\^2i " ^\i) Jt 

It follows from the similar argument of fl2.12p that 

dF\ ^ c^XijMi 1 ^ c^AaiM, 1 ^ 2c^M, 

c^t di\\2i — Aij) Aij di{X2i — Aii) A2i di{X2i — Aij) 

It follows from flTT^ that 

Therefore |^(y9, y9^)(t)|p < ^ .^"^ j|^^^ . for alH = 1, ■ ■ ■ ,n. Hence F is equicon- 
tinuous on C[o,k](IR, I^")- In view of Lemma [2. 2 [ F(r) is uniformly bounded. 

Next we claim that F : F — )■ F is compact. Define the operator sequence 
{F(")}, where F^") : C[o,k](M,M") -^ C[o,K](ffi,K") by 

{F{^,^r){-n), te{-oo,-n), 
F(^,^,)(t), te[-n,n], 
F{ip,ipr){n), t e (n, +00). 

Hence the sequence {F^")} are uniformly bounded and equicontinuous. It follows 
from Arzela-Ascoli theorem that F*^"^ is compact. Therefore we have 

|f/") - F4 = sup \Fp{^,^rm - F,{^,v,){t)\e-^\'\ 

sup |i^"(^,^.)(t)-F,(yp,yp.)(t)|e-^l*l (2.17) 

tG(— 00,— n)U(n,+oo) 

< 2Kie~f^ ^ as t ^ 00, 



where ^9 G C[o,k](1R5 1R")- By the virtue of proposition 2.1 in [32], the sequence 
{F'-"^} converges to F in F with respect to the norm | ■ \p. Therefore F : F — )• F 
is compact. 

D 

Theorem 2.1. Assume that (Hi) and (H2) hold. Suppose that(p E C[o,ii:\(M.,W^) , 

(p G C[o,K](1^5l^'^) ^e a pair of upper and lower solutions of h2.2) . and 

hm^(t)=0, hm^(t) = K, (2.18) 

t— >— 00 t— >+oo 

then li2.S\) and ^2.4^ admit a solution. That is, the problem Ii2. 1\) has a traveling 
wave solution. 

Proof. First we define the following profile set such as in Lemma 12.21 

r = {v^GC[o,K](K,K"): 0<^<^}, 
it is easy to show that F is a closed convex set. 



Now we define the operator such as in Lemma 12.31 

F : F ^ F, 

in view of Lemma 12.21 Lemma 12. 3[ Lemma 12.41 F is continuous and compact 
with respect to | ■ \p. By Schauder's fixed point theorem, there exists a fixed 
point (f* eT such that F{(f*, (f*) = (f*. Using the Perron Theorem, (12. 8p has a 
solution (f*, that is (f* is a solution of (12. 2p . 

Next we will show that (f* satisfies the boundary condition (12.41) . In view of 
Lemma 12.21 (p < (p* <'^. It follows from (I2.18P that 



< lim (p*{t) < lim ^ = 0, 
t— >— 00 t— ^— 00 

K > lim if*{t) > lim (^ = K. 

Therefore ip* is a traveling wave solution of the problem (12. ip . Thus the proof 
is completed. D 

3 Coupled quasi-upper and quasi-lower solu- 
tions 

In Theorem 2.1, we see that the smooth conditions on the coupled upper and 
lower solutions are too strong. In fact, it is very difficult to seek the C^ smooth 
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coupled upper and lower solutions for a special model. We intend to relax the 
smoothness of the upper and lower solutions to C^. Thus we should cite the 
modified Perron theorem in |2]. 

Definition 3.1. Considering the following scalar ordinary equation 

u"{t) + au{t) + Pu{t) + f{t) = 0, teR, u{t) e M, (3.1) 

where P < 0, f is a bounded and continuous function onR\ {0} and both f{0~^) 
and /(0~) exist. Then, a function u defined on M is said to be a generalized 
solution of h3.1\) if 

(1) u and u' are bounded and continuous on M. 

(2) u" exists and is continuous on ]R\{0}, and both n"(0^) andu"{0^) exist. 

Lemma 3.1. (i^) Consider Ii3.1\) with P < 0, and assume that 

(1) f is a bounded and continuous function onM.\ {0} and both /(0+) and 
f{0^) exist, 

(2) 113. 1\) holds in the classical sense for all t except possibly at t = 0. 
Then liS. 1\) has a unique generalized solution u given by 

1 rt r+(x> 

u{t) = ( / e^'^'^'^f{s)ds+ e^'^'-''>f{s)ds), (3.2) 

where Ai and X2 are respectively the negative and positive roots of X^+aX+P = 0. 

Now we give the following definition of an admissible upper and lower solu- 
tions. 

Definition 3.2. Assume that ^, (p e Cl{R,W) , §(t) and ^(t) exist and 
continuous onR\ {0}, and 

,d^Lp , ,, ,. d^(p . , ,. d^Lf) , , 

sup — --(t) < +00, and lim — — (t), lim — — (t) exist, 

t^R\{o} dt^ t^o- dt^ t^o+ dt^ 

.d^(p , ,, ,. d'^ip , , ,. d^(p , , 

sup K-7(t) < +00, and lim — -^-(t), lim — -r-(t) exist, 
t^M.\{o} dt^ t^o- dt^ t^o+ dt^ 

(p, Cp satisfy 

m)-§'f'Iit)>M<fi, ma,, [0k, [VrU, [0rU),for all teR\{0}, 
c0'i{t) - %0'm < f^{0,, [pU, [^]6„ [prU [0rU),for all teM\{0}.^ " ^ 

for i = 1,- ■ ■ ,n. Then ip and are called coupled quasi-upper solution and 
quasi-lower solution of Ii2.2\) . respectively. 
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Lemma 3.2. Ifip, ^ C[o,k](I^5I^") O'I"^ o, pair of quasi-upper and quasi-lower 
solutions of \2.S\) . Then x*-^^ and x*-"^-* constructed by \2.1C^] . where (p and (p are 
replaced with these quasi-upper and quasi-lower solutions, are a pair of upper 
and lower solutions of ^2.2\) . moreover 

< X^^) < x(^) < If. 



Proof Combining (13 ■4p and fl2.10p yields 

i^\^^y - -iix[^^)" + A^r^ = Pm + M'fi, ['fU, [^]fe,:> [VrU, [^r]dj 



C2" 



<^'i- 4^i + Pi^u for all te^\ {0}. 



In view of Lemma 13.1^ we have 

^''(t) < ^ ( f e^"(*-)(^:(.) - ^^r(.) + I3m{s))ds 

+ n e''-^'-^\0[is) - %0'l{s) + ^0.{s))ds) ^^-^^ 

Jt c 

= (pi{t), for all t e R\{0}. 

In a similar way, we have x^^ (t) > 0i{t) for all t G M \ {0}. Since that x*-^-*, 
x(^) satisfy (l^3U]) . x^^), x^^) e C[o,k](K,M") n C^2(m^M") and ([13]) holds for all 
t G M. A similar argument in Lemma 2.1 of [22| shows that x*-^-*, x*-^-* are a pair 
of coupled upper and lower solutions of (12. 2p . D 

Theorem 3.1. Assume that {Hi) and {H2) hold. Suppose that (p G C[o,k](I^; M"), 

(p G C[o,K](I^5l^") be a pair of coupled quasi-upper and quasi-lower solutions of 
HJj, and 

lim (p{t) = 0, lim (p{t) = K. 

t— >— 00 t— >+oo 

Then 1(2. ^)) and ^2.4^ admit a solution. That is, the problem 112. 1\} has a traveling 
wave solution. 
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Proof. Let 



^2 



^^ = 77X r^( / e "^ nA^. + M^^, [^]a„ [^]6., [^rU, [^rU))ds 



"+00 



+ / e^'^'^'-'\l3i^i + fi{^i,[(p\aA0\h, [<^r]c„ [<^r]dj)0?s), 



^^^'^ = -7n:^^-\^ f e^^^^'-^\(5,(p, + m,, [^],,, [^],„ [^.],,, [^.],J)rfs 



—00 
00 



for i = 1,2, ■ ■ ■ ,n. Then by Lemma 13.2^ x'-^^ and x'-^-* are a pair of coupled 
upper and lower solutions of (12.21) . Replace ip, (p with x'^^^ x^^\ We then use 
Theorem 2.1 and obtain the same results directly. D 

4 Ordered quasi-upper and quasi-lower solu- 
tions 

If f ((/3, (pr) is quasimonotone nondecreasing, that is, 6j = dj = for all z, then 
the upper and lower solution are ordered. In this case [21 ESI ES] have showed 
the existence of the traveling wave solution. 

Definition 4.1. A pair of vectors (p = ((^i,--- ,(pn)i = {0i,''' ,0n) in 
Cl{M.,W^) are called ordered upper and lower solutions of 1^2. S\} if (p > (p and if 

0m - 7i0Kt) < M0, 0r), {^ = l,■■■,n), ^ ' ' 



where iprit) = ip{t — r). 



Definition 4.2. Assume that (p, e Cl{W.,W'') , §(t) and ^(t) exist and 



continuous on'K\ {0}, and 

sup — --(t) < +00, and lim — — (t), lim — — (t) extst, 

f_j.iR\{o} dt^ t^o- dt^ t^o+ at^ 

,d^(p d'^ip , , ,. d'^ip . , 

sup — --(t) < +CX), and lim — — (t), lim — — (t) extst, 
i^iR\{o} dt^ t^o- dt^ t^o+ dt"^ 

(^, (f satisfy 

<fi{t)-^0i{t)>M0, Vr),for all teR\{0}, 
'f'iit) - ^0i{t) < M0, 0r)jor all tGM\{0}, 



(4.2) 



(4.3) 
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for i = 1, ■ ■ ■ ,ra. Then ip and are called ordered quasi-upper solution and 
quasi-lower solution of h2.2^) . respectively. 

Next theorem shows that if the upper or quasi- upper solution (p{t) is non- 
decreasing respect to t, then the solution ^*{t) is also nondecreasing. Using 
the similar argument in Lemma 4.1 of [22], we can induce the existence of the 
traveling wavefront of fl2.2p . 

Theorem 4.1. Assume that {Hi) and {H2) hold andf{ip,ipr) is quasimonotone 
nondecreasing. Suppose that (p G C^o.kjIIRjIR"); V G C[o,k](I^5I^") ^e a pair of 
ordered quasi-upper and quasi-lower solutions of /i2.^) . and (pit) is nondecreasing 
with respect to t, 

lim ip{t) = 0, lim 0{t) = K. (4.4) 

Then the solution of l\2. ^) and ^2.4^ ip* is nondecreasing with respect to t. That 
is, the problem 112. 1\} at least has a traveling wavefront solution. 



In Theorem 4.1, the condition (14. 4p can be replaced by further restrictions 
on f , that is 

(HI) f(u, u^)|u=o = f(u, u^)|u=K = and f(u, u^)|u=l 7^ for the 

constant-valued function L with < L < K and L 7^ 0,L 7^ K. 

As similar as the argument in Theorem 4.4 of [22], Theorem 4.1 can be 
transformed into 

Theorem 4.2. Assume that {H^) and {H2) hold and i{ip, ipr) is quasimonotone 
nondecreasing. Suppose that ip G C[o,k](IR, IR"); ^ C[o,k](I^!I^") b^ a pair of 
ordered quasi-upper and quasi-lower solutions of Ii2.2\) . and (p{t) is nondecreasing 
with respect to t, 

< 0{t) < <p{t) < K, (p{t) ^ 0, !p{t) ^ K in M, (4.5) 



then the solution of l{2.2\) and lj^2.4\) ^* is nondecreasing with respect to t. That 
is, the problem Ii2. 1\) at least has a traveling wavefront solution. 

Remark 4.1. In Theorems 4-i o-nd 4-2, the ordered quasi-lower solution is not 
necessary nondecreasing. The results of Theorem 4-1 has been obtained by many 
authors, for example. Theorem 2.2 of 12^ and Theorem 11 of 12], and Theorem 
3.6oflJB]. 
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In the similar way, Theorem 3.1 can be transformed into 

Theorem 4.3. Assume that (H*) and {H2) hold. Suppose that(p e C[o,k](IR, ^^), 
(p G C[o,K](I^;I^") ^e a pair of coupled quasi-upper and quasi-lower solutions of 

< ^(t) < (p{t) < K, ^(t) 7^ 0, (p{t) ^ K in M, (4.6) 

Then /i2. ^) and ^2.4\ ) admit a solution. That is, the problem ^2. 1\) has a traveling 
wave solution. 



5 Application 

In this section, we use our results obtained in previous sections to consider the 
delayed reaction diffusion models. 

5.1 Belousov-Zhabotinskii equations 

Consider the delayed Belousov-Zhabotinskii equations 



^^^ - ^^S^ = u{x,t){l-u{x,t) -rv{x,t-T2)), 
^^^ - 4S^ = -bu{x, t - Ti)v{x, t), 



(5.1) 



where u{x, t),v{x, t) are scalar functions, r > 0, 6 > are constants. In the bio- 
logical sense, u and v represent the Bromic acid and bromide ion concentrations 
respectively (see more details in J27j). Without the delays, the existences of the 
traveling wave solutions were considered in [HI [T21 E31 ES] • When ri = 0, r2 7^ 0, 
[21 |23l [36] studied the traveling wave solution by using of various methods to 
construct quasi-upper solution. 

Now we seek the traveling wave solution, whose form is u{x, t) = (fi{t-\-x/c), 
v{x,t) = ip2{t + x/c). Clearly the wave equations corresponding to (12. 2 p is the 
following form 

^[(t) - ^ip'lit) = ^i(t)(l - ^i(t) - rMt - ^2)), ., ^. 

^m - ^^m = ~bMt)Mt - ri). ^''■^' 

We seek a traveling wave solution of (15. ip with the boundary conditions 

lim flit) = 0, lim f2(t) = 0, 
t— >— 00 t— >— 00 

lim fi{t) = 1, lim if2{t) = 1- 

i— >'+oo i— >+oo 
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It is easy to check that {H^) and {H2) are satisfied. We only need seek the 
coupled quasi-upper and quasi- lower solutions of f l5.2p . 

li {p = (^1^1,(^2) and if = ((^1,(^2) are coupled quasi-upper and quasi-lower 
solutions of fl5.2p . they must satisfy 



^[{t) - ^^'/(t) > (^i(t)(l - (^i(t) - r02{t - T2)), for all t eR\ {0}, 
^2it) - Wiit) > -h(p2{t)^i{t - Ti), for all t G M \ {0}, 
^[{t) - ^0'l{t) < (^i(t)(l - 0i{t) - r^2{t - T2)), for all t G M \ {0}, 
(p'2{t) - ^if'iit) < -b02{t)^iit - n), for all t G R \ {0}. 



(5.3) 



Assume that ip and (p 

^.(*) = { f_ 1,-M. ; ; „: -^ ^^(*) = { f_ I'e-., t > 0, (") 

where Ai, A2, A3, 5, fc are undetermined positive constants. 
Direct calculations show that 

^'■W = { |:-.\ *;o: and ^UO = { 1:-.. l-ll (5.6) 

^^ \, ' and ^^'(t) = <^ 2^ ' ^-^' (5.7) 

., ,.. / A3<5A;e^^*, t < 0, , .„.,. / Ai^fce^^**, t < 0, ._ „. 

^^^'^ = [ Xs5ke-'^\ t > 0, ""' ^i(') = 1 -Xl6ke~'-^\ t > 0, ^^'^^ 

From 05.6p .( 13T|) and (15. 8p . we see that the first derivatives are continuous and 
the second derivatives exist and continuous on M \ {0}. Hence ip and satisfy 
(14. 2p . Now we will choose proper Ai, A2, A3, 6, k such that (15. 3p holds. 

Substituting (15. 4p . (15. 6p and (15.71) into the first equation of (15. 3p yields 



2 c2 2 -2^2 

-e-"^* - i(-^e-^^*) > (1 - -e- ,,. . , 
2 c2^ 2 ^ - ^ 2 '' 2 



2 ' ; 1 (5-9) 

e"""'* - ^(-— e"^'*) > (1 - -e"^'*)(l - 1 + -e^^^*) for t > 0. 
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In order to induce fl5.9p . the following is sufficient 



\2 

Ai - ^ - 1 > 0. (5.10) 



Hence we set 



cH\-J\-^ 



r-2. 



Ai = \ . (5.11) 

Substituting (I5.4p -f l5l8|) into the second equation of (15. 3p yields 

^e^^* - \%^-' > -b^e'^'Ske'^^'-^^^ for t < 0, 
2 c2 2 - 2 

^e-"^* - 4(--e"''*) ^ -^(1 - -e-^''')6ke^-'^'-^^^ f or < i < n, (5.12) 

Zi i^ Zi Zi 

^e-"^* - ^(-^e-^^*) > -6(1 - \-^-%k - 5A;e-"«(*-^^)) for t > r,. 
2 c2^ 2 ^ ~ ^ 2 ^^ ^ 

In order to induce (I5.12p . the following is sufficient 

A2-4>0- (5-13) 

Hence we set 

A2 = ei, where £1 << 1. (5.14) 

Substituting (I5.4p -( 13T5]) into the third equation of (15. 3 p yields 

XsSke^^' - \xl6ke^'' < 6ke^''\l - 6ke^'' - r-e^^^*-"^)) for t < 0, 

XsSke-^^' + \xl6ke-^^' < {k - 5ke-^-''){l -k + 6ke-^^' - r^e^^^*-"^)) 
c^ 2 

for < t < r2, (^-1^) 

\2 1 

4) < (A; - (5A;e-^«*)(l - fc + Ske-^'' - r(l - -( 
c^ 2 

for t > T2. 



Ske-^''\X3 + -§)<{k- 6ke~^''){l -k + 6ke-^^' - r(l - -e"^^*)) 



If we set 6 « 1, A2 << 1, in order to induce (I5.15p . the following is sufficient 



\2 

A3 - ^ - K 0, 



5ke-^'\X; + ^ + l-k)<k{l-k), (5.16) 



c2 



\2 

Ske-^'\X3 + ^ + 1 - r ~ k) < k{l - r - k). 



c2 
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To ensure (15.1 61) . the parameter must satisfy 

r < 1 (5.17) 

and we set 



c^(i-Vi-^) 

As = ^ £2, where £2 << 1,/i; << 1- (5.18) 

It is easy to see that the fourth equation of (15. 3p naturally holds. Therefore, we 
have proved 

Lemma 5.1. Assume that the parameter of the problem Ii5.1\) r < 1. Then 
there exists a constant c* = 2 such that if c> c* , ip and (p, which are defined in 
( [5.^[ ) and \5.J^ , are a pair of coupled quasi-upper and quasi-lower solutions of 
the problem l\5.^) . 



Finally, by Theorem 4.3 we have 

Theorem 5.1. Assume that the parameter of the problem 115. 1\) r < 1. Then 

there exists a constant c* = 2 such that if c > c* , the problem Ii5.1\) has a 

traveling wave solution u{x,t) = ipi{t + x/c), v{x,t) = {p2(t -\- x/c), which 
connects (0,0) and (1,1). 

Remark 5.1. In J^, the critical value of wave velocity c* is dependent on the 
parameter b. In our Theorem 5.1, we show that c* may be a constant. 

5.2 Mutualistic Lotka-Volterra model 

The delayed mutualistic Lotka-Volterramodel is as follows 



^^^ - di ^ f^/^ = ru{x,t){l - aiu{x,t) + biv{x,t)), 
^^^ - d2^-^^ = v{x, t){a2u{x, t - r) - 62), 



(5.19) 



where u{x,t), v{x,t) are scalar functions, and r, ai, 02, &i, &2 are all positive 
constants, di,d2 are the positive diffusion coefficients, r represent the positive 
delay. For a detailed description of this model, we refer to [27]. The above 
model with r = 0, 61 < has been considered in [U O [27]. The case with 
r = 0, 61 < was studied in [21 [23]. However, when 61 < the reaction term of 
(I5.19P is not quasimonotone nondecreasing, there was not sufficient condition 
to prove Lemma 3.5 of [23j and Lemma 12 of [2]. 



If we assume that 

a2 > ai62, (5.20) 

then the model (I5.19P has a unique positive equihbrium 

K,--) = (^.^(^-l)). (5.21) 

^2 Ol ^2 

The wave equations corresponding to (12.21) is the following form 



(5.22) 



ip[{t) - ^ifi'lit) = r^i(t)(l - aiy^i(t) + hM^)), 
V'2it) - ^^'iit) = Mt)ia2Mt -r)- 62). 

We will find the solution of the above wave equation such that 

lim v^i(t) = 0, lim 'f2(t) = 0, 
t— >— 00 t— >— 00 

lim (pi{t) = u*, lim (p2{t) = v*. 



It is also easy to verify that {H^) and {H2) are satisfied. We only need seek 
the coupled quasi-upper and quasi- lower solutions of (15.221) . 

li if = ((^1,(^2) and = ((^1,(^2) are coupled quasi-upper and quasi-lower 
solutions of (I5.22p . they must satisfy 

if[{t) - ^ip'lit) > r^i(t)(l - aiifiit) + h02{t)), for all t G R \ {0}, 
if'^it) - ^0l{t) > <^2(t)(a2^i(t - r) - 62), for all t G M \ {0}, 
'P'lit) - ^'P'lit) < r0i{t){l - ai0i{t) - bi02{t)), for all i G M \ {0}, ^ ' ^ 
02it) - %02{t) < 02{t){a20i{t - r) - 62), for all t e^\ {0}. 

Assume that (p and 

( HlpAii f < n f Hlp^2i f < 



^'W = { f-'^L->'. t>0: -"^'^W-O. (5.25) 

where Ai, A2, A3, 5, k are undetermined positive constants. 

As similar as the process in the above subsection, in order to (I5.23p . the 
following is sufficient 
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Ai- 


c2 


-r > 0, 


Ai + 


diXl 
c2 


- hiv* > 


A2(l 


- %\.) > 0, 


As- 


M. 


r, 


k «l,5 ■ 


« 1. 



(5.26) 



If 



c>max{2Vrrfi,A/rfi(r + ^^-l)}, (5.27) 

0-2 



set 



cm + Jl-^r^: 



Ai = ^^- , As « 1, A3 « 1, k«l,6 « 1, (5.28) 

then the sufficient conditions (I5.26|) hold. Therefore by Theorem 4.3, we have 

Theorem 5.2. Assume that the parameter of the problem i\5.19\) satisfy ^5.20\} . 
Then there exists a constant c* = lalax{2^/rd^, \/di{r + ^^ — 1)} such that if 
c > €*, the problem /i5. 1]) has a traveling wave solution u{x,t) = ipi(t + x/c), 

' hi 1 l a\h2 
^ 02 ' 61 ^ a2 



v{x,t) = ^2{t + x/c), which connects (0,0) and (^, f (^i^ - 1)). 



6 Discussion 

We aim to study the existence for the travehng wave solution of the discrete- 
delayed reaction diffusion systems, where the reaction term is mixed quasi- 
monotone. Our result is that the existence of the coupled quasi-upper and 
quasi-lower solutions ensure the traveling wave solution exists. In the equa- 
tions of [21 |23l |36], the reaction term is quasimonotone nondecreasing. In fact 
the conditions of quasimonotone nondecreasing property is very strong. The 
predator-prey model which is studied in [21 123] do not satisfy the quasimono- 
tone nondecreasing property, and the existence theorem of [21123] is not suitable. 
Hence the application scope of Theorem 3.1 is wider than [21 [231 [36]. Moreover, 
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as a special case of Theorem 3.1, Theorem 4.1 contains the previous existence 
theorem for the travehng wave solution. 

Our technique to deal with the mixed quasimonotonity is constructing the 
coupled upper and lower solutions. Recently the method so-called cross iteration 
scheme was developed in p!9] to deal with the traveling wave solution for the 2 
dimensional competitive Lotka-Volterra model. Comparing with the model of 
[T9] . the systems in this paper are extensive to n dimension. 

The classical coupled upper and lower solutions need the second order smooth- 
ness. It is very difficult to satisfy this condition in the real model. Thus it is 
necessary to relax the smoothness of the coupled upper and lower solutions 
to first order smoothness, which is called coupled quasi-upper and quasi-lower 
solutions. This paper apply the modified Perron theorem with the case C^ 
smoothness, which is first proposed in [2] . Our existence theorem of the travel- 
ing wave solution is suitable to all 2 species or 3 species Lotka-Volterra systems. 
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